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Abstract 


Multi-scale problems in numerical electromagnetics are becoming increasingly common with the advent and widespread usage 
of compact mobile phones, body area networks, small and nano antennas, sensors, high-speed interconnects, integrated circuits 
and complex electronic packaging structures, to name just a few commercial applications. Numerical electromagnetic modeling 
and simulation of structures with multi-scale features is highly challenging due to the fact that electrically small as well as large 
features are simultaneously present in the model which demands for discretization of the computational domain such that the 
number of degrees of freedom is very large, thus, levying a heavy burden on computational resources. The multi-scale nature of 
a given problem also exacerbates the challenge of generating very fine meshes which do not introduce instabilities or ill- 
conditioned behaviors. 


In this work we introduce a hybrid technique, which combines frequency domain and time domain techniques in a manner such 
that the fine features (electrically small) of the object being modeled are handled by the Method of Moments (MoM) technique 
while the electrically large parts of the structure are dealt with by using the Finite-Difference Time-Domain (FDTD) technique in 
order to reduce the computational burden. Recently, structures with multi-scale features have been simulated by using the dipole 
moment (DM) approach combined with the FDTD technique to handle fine features in a multi-scale geometry. However, when 
the size of the scatterer becomes larger in terms of the wavelength and the quasi-static assumption becomes invalid, extensive 
modifications of the DM/FDTD hybrid approach are needed resulting in a high computational cost. 


The research proposes a novel hybrid FDTD technique, which combines the Method of Moments and the Finite-Difference Time- 
Domain techniques directly in the time domain circumventing the need to carry out frequency transform calculations as required 
in the DM approach when the object size is not small (size>A/20). The proposed technique utilizes piecewise sinusoidal basis 
functions to represent the currents on arbitrarily shaped wires with fine features, and modified RWG basis function for surfaces. 
The fields scattered by the object with fine features in MoM region are computed in the time domain on a planar interface. The 
time domain fields obtained at the planar interface are then combined with the FDTD update equations. In contrast to the 
existing techniques used to handle this type of problems, the proposed technique is both efficient as well as stable. 


Index Terms: MoM, FDTD, DM, RWG, DoF, TDIE, EFIE, MFIE, CFIE, Matrix, RHS, Finite Difference, Lagrange Quadratic 
Interpolation, Modified RWG 
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Presentation Outline 


Introduction 


e Research motivation 


Multi-scale problem definition 





Integrated circuits 


antennas 


Novel techniques to handle multi-scale problems 


* Straight wire problem 





External 
antenna 


rend 


e Arbitrary shaped wire problem 





Internal 
antenna 





e Wire and surface problem 





FDTD region 


MoM region 


* Capacitance matrix computation for coated multi-conductors 


Observations and Conclusions 
e Future research 


: References 


42 cells 





V=0 boundary condition is used 
on all four sides 
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Introduction: Modeling Challenges 1n Current 
Technologies 


e Ubiquity of multi-scale problems ın numerical electromagnetics. 


Conformal 
printed 
antennas and 
circuits 







———— 4 





Body area 
networks 


MRI scanning using rapid 7T coil 
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Research Motivation 


* Simulation of multi-scale problems with fine features is highly 
challenging. ; ee eae 


Scenario 1: Utilization of FDTD method 
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Scenario 3: Utilization of FEM and MoM method 
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Amplitude Comparison of E., from New Conformal FDTD 


Scenario 2: Utilization of FEM method et -——— New Conformal EDTD 
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Fine feature of object captured by à s as ae 
s g^ TR U^ 
conventional FDTD mesh. Extremely fine mesh results in high E S 

computational cost. an _— 

n "n — a E Utilization of fine mesh results in 0.6 d 
ghest dielectric constant. instability and does not guarantee 
accurate results using different 0 02 04 06 08 1 12 14 16 18 2 


type of solvers. Distance along Y in 


E 
Novel hybrid technique to handle S 


multi-scale problems 





Time-domain scattered fields at the Fine features are handled by the 
interface are used as a source for the , . I 
FDTD domain time-domain compatible MoM code 


FDTD domain ——— Planar Interface 


Mesh size: 4/20 MoM domain Time-domain compatible basis functions are used 


i @ highest frequency of interest to represent current on geometry with fine 


features. These basis functions are used to obtain 
= Antenna with fine , , ; . 
faures ia oM doman scattered field at an observation point directly in 


Object with electrically large 
features in FDTD domain 





Soupes dade Tossantemis the time-domain without using Inverse Fourier 


MoM domain are obtained Transform 
directly in time-domain. 
The novel hybrid FDTD technique combines 


the MoM and FDTD techniques directly in 
the time domain. Time-domain scattered fields are sampled at a 





Efficiency factor for the hybrid method is large compared to the planar interface 


previous approach where scattered fields at interface nodes were 
computed in frequency domain, and inverse transform was used to 
obtain fields in the time-domain. 


Sampled time-domain scattered fields are then used as 


source for the FDTD domain 





Time step of the time-domain compatible MoM technique is 
chosen such that it is equal to the time-step of the FDTD technique. 


Time domain compatible MoM code is different than 


, , l l the TDIE formulation, and is stable. 
Different types of problems investigated are discussed in the 


following sections. 


Straight wires 
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Time-domain compatible basis functions 
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Frequency-domain scattered field expressions: 
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Time-domain scattered fields at the 
interface are used as a source for the 
FDTD domain 


ME Planar Interface 
MoM domain 
\ "—. Antenna with fine 
features in MoM domain 


Scattered fields from antenna m 
MoM domain are obtained 
directly in time-domain. 


FDTD domain 


Mesh size: 3/20 


À (à) highest frequency of uA 


Object with electrically we. 
features in FDTD domain 





Time-domain compatible MoM formulation is different 
as compared to the conventional TDIE formulation, and 
is stable. 


The basis function enables us to obtain scattered field 
directly in the time-domain at interface nodes without 
using the Inverse Fourier Transform, this is in contrast 
to the previous approach. 


The number of operations required to obtain time- 
domain fields at interface nodes is significantly reduced 
compared to past approach. 


cospoH approximation in the frequency domain makes 
scattered field expressions time-domain compatible by 
avoiding causality violation in the time-domain scattered 
field expressions. 





NUMERICAL RESULTS 


Hybrid FDTD method for problems involving 
straight wire antennas 
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Validation of ‘Time-Domain compatible 
MoM Code 


Note: The wire can be tilted or have 
an arbitrary curved shape 


X 
à @ highest frequency of interest: 4 GHz 








In the frequency domain, the Gaussian 
pulse has a range of 50 MHz to 4 GHz 
with 50 MHz frequency step 


a 


Straight wire antenna of length 0.75 cm with 
a delta gap voltage source 


Observation Line (Cosine modulated Gaussian pulse) 
A straight wire is directed along x. 


Length of the straight wire is 0.75 cm. Wire antenna radius — 0.005*length 


The radiated fields are computed in the time-domain along the observation line, and transformed 
to the frequency domain, so as to compare results at a desired frequency with commercial MoM 
solver to show results in time-domain using time-compatible MoM code are accurate. 


Results along the observation line are compared at 2 GHz. 


PENNSTATE 


Time-Domain compatible MoM Code vs FEKO 
at 2 GHz 
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In the frequency domain, the Gaussian 
pulse has a range of 50 MHz to 4 GHz 
with 50 MHz frequency step 


2 @ highest frequency of interest: 4 GHz 








Z 
F- Straight wire antenna of length 0.75 cm with 
a delta gap voltage source 


Observation Line (Cosine modulated Gaussian pulse) 


A straight wire is directed along x. Wire antenna radius — 0.005*length 
Length of the straight wire is 0.75 cm 


The radiated fields are computed in the time-domain along the observation line, and transformed to the frequency domain. 
Results along the observation line are compared at 2 GHz which is the center frequency of the modulated Gaussian pulse. 
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Hybrid FDTD Technique 
Example: Tilted Wire Antenna in Transmit Mode 


; FDTD domain mesh size: 2/20 
Variable Value | 
A @ highest frequency of interest: 8 GHz 










f The modulated Gaussian pulse spectrum 
À 3. /9 CM ho @ center frequency : 4 GHz has a range of 40 MHz to 8 GHz with 40 
MHz frequency step. 
À /.5 cm 
0 n 


Interface 


A/20 0.1875 cm (26.25 cm x 26.25 cm) 


MoM domain 
Tilt angle=45° 
Observation Line 


FDTD domain 
| 


e——_—_—__—— IL 






X 
f Straight wire antenna of length 0.1875 cm 
with a delta gap voltage source 


(Cosine modulated Gaussian pulse) 
0.1875 cm 


Wire antenna radius - 0.005*length 


pec plate 
pec plate touches FDTD domain 
pml boundaries along x and y 
directions which renders it infinite 





MoM domain radiated fields at interface 
^ — nodes are obtained in the time-domain 


Cross-sectional view 


A straight wire antenna is inclined to x axis at an angle=45° in xz plane. 


Length of the straight wire antenna is 0.1875 cm. 
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E, Comparison at 4 GHz 
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Results for radiated fields along 
observation line are in close agreement 
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Effici factor = 
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Croce sectional view 
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Distance Along Z in A 


FDTD domain mesh size: 3/2 
à. i$ highest frequency of interest: 8 GHz 


io € center [requency: 4 GHz Tbe modulsicd Gaussian pulse specinam 
= 7 X has a imge ofdo MHz lo 8 GHz with 40 
x MHz 







Interface 
(26.25cm x 26.25 em) 


Moh demain 
Tilt angle=45° 


bt; 


=. ^ 
= Sirake wire amber of lcmguh 0.1375 cm 
with adele gap voltage wource 


i 4 GHz is the center 

ne (Cosine medelwied Cured pulse) 

— = ¢,606 "length frequency of the 
menm o modulated Gaussian 

pulse 


Observation Line 
4X 


pec piste Touches FTETTE domain 
pad boundane: along x andy 
directions which renders it intiwibe 


À straight wire antenna is inclined to x axis at an angle-45" in xz plane. 


Length of the straight wire antenna is 0.1875 cm. 








Hybrid FDTD Technique T 
Example: Scattering from a Tilted Wire 


à @ highest frequency of interest: 8 GHz 


À S cm ho @ center frequency : 4 GHz In the frequency domain, the modulated 


Gaussian pulse has a range of 40 MHz to 
Ào aoe ii 
| FDID domain 
d 
| 7.5 cm 








8 GHz with 40 MHz frequency step. 
)/20 0.1875 cm MoM domain 


Interface Tilt angle=45° 
(26.25 cm x 26.25 cm) 





Observation Line k í E . 
Incident modulated 


Z gaussian pulse 





^x 





R Straight wire antenna of length 0.1875 cm 
| 
| 


lat 
pee prawe 0.1875 cm Wire antenna radius - 0.005*length 
\ MoM domain scattered field solution at 
—>Y pec plate touches FDTD domain pml interface nodes is obtained in the time-domain 


boundaries along x and y directions 


Cross-sectional view which renders it infinite 


A straight wire is inclined to x axis at an angle=45° in xz plane. 


Length of the straight wire is 0.1875 cm. 
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E, Comparison at 4 GHz 


Amplitude of E. Phase of E. 
200 ' | 









—*— |mage theory 
—— Hybrid FDTD 


5 8 


Amplitude in V/m 
Phase in Degree 
on 
eC 





“2 4.5 4 -0.5 0 


Distance Along Z in A Distance Along Z in A 


Mesh size: 4/20 


4. @ highest frequency of interest: 8 GHz 


ig (@ center frequency: 4 GHz ee a Deo 
Gaussian pulse has a range of 40 MEE to 
X & GHz willi 40 MHz frequency step. 
| MoM domain 
FDTDdomain | Interface . Tilt angle=45° 
| (26,25 cm x 26,25 cm) i 


one «J -i Incident modulated 
Scattered field results along d = — i puce pert 
i wire antenna of le 0,1575 cm . 
Dopo v eges 4 GHz is the center 


0.1575 cm Wire antenna radius = 0.005*length 





observation line are in close agreement 
frequency of the 





MoM! domain scattered field solution is 


—»V pec plate touches domain pm obtained in the time-domain modulated Gaussian 
boundaries along x and y directions 
Number of operations using LE.T. method at interface 7 Crosssectional'view which renders it infinite 
10 A straight wire is inclined to x axis at an angle=45° in xz plane. pulse 


Efficiency factor = F - F 
i Number of operations using hybrid method 
Length of the straight wire is 0.1875 cm. 


PENNSTATE 





Handling arbitrary shaped wires with bends 


Scattered field expressions from basis function 


Branch 1 current ü 
current: 


^ distribution along u, / 
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Bent wire and basis function current on the branches 
Virtual current approach 
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NUMERICAL RESULTS 


Hybrid FDTD method for problems involving 
arbitrary shaped wire antennas 


Problem definition “oe 
Validation of the MoM code 





Tilted Circular Loop 


In the frequency domain, the modulated 
Gaussian pulse has a range of 40 MHz to 


i 8 GHz with 40 MHz frequency step. 
Interface tire: GA, t GA. 02272 B cm 225 cm 


| Plane wave 
Noder on interface have a tpacing of A... (20 3 
in yands directions, Theta =45, phi=45 deg. 
121 x 121 nodes on tha Interface. 


Mq OB GH: [highest irequancy al interest) 






Loop is centered 

at origin (0,0,0) 

and tilted by 45° 
about x axis 

Loop diameter: 

3 fdtd cells(5.625 mm) 


Interface scattered fields obtained using iu 
modified DM approach and FEKO. ( ite iinet 


interface is parallel to yz plane at a distance of si | n corresponds to z-0 
0.75 cmfrom the center of loop along x direction. Incident modulated Loop diameter: 5.625 mm. 





1 
gaussian pulse Loop cross section radius = S08 of the loop circumference. 


Scattered fields are computed in the time-domain along the 
observation lines, and transformed to the frequency domain. 


MoM code vs FEKO results comparison E 
Observation line | 


Obs. Line 1: Interface E, field comparison Obs. Line 1: Interface E, field comparison 





Magnitude Phase Magnitude Phase 


Ey magnitude (vim ) 
Ey phase (disgree) 





Interface Nada 





Line 1: x = 0.75 cm, z = 0 cm 
Tilted Circular Loop 


In the frequency domain, the modulated 
Gaussian pulse has a range of 40 MHz to 


leterfuco sine GA, «OL 225 ome 225 om 8 GHz with 40 MHz frequency step. 







e^ inerfuce heve a pacing of Ap /20 3 Plane weve . 
in y vd t Grections, Theta = 45, phi=45 deg. 
221 & 122 node: on the interface. 
din @8 GH (highent frequency of interest) 2 
” m 
p 
Loop is centered Fs 
at origin (0,0,0) 
and tilted by 45° 
about x axis 
3 fdtd celis5.625 mm) 2 
interface scattered fields obtained using 
modified OM approach and FEKO. ( Obs. Une 1 
interface is parallel to yz plane at a distance of 4- corresponds to z=0 
0.75 cm from the center of loop along x direction. Incident modulated a ee 
gaussian pulse Loop cross section radius © of the loop circumference. 


Results along the observation line are compared at 4 GHz which is the center 
frequency of the modulated Gaussian pulse. MoM results agree well with FEKO results. 
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MoM code vs FEKO results comparison 
Observation line 2 


Obs. Line 2: Interface E, field comparison Obs. Line 2: Interface E, field comparison 


Phasa 


Magnitude 


Phase 


Magnitude 














Line 2: x = 0.75 cm, z = -3.75 cm 


Tilted Circular Loop mm 
Gaussian D a range of 40 MHz to 
Itecfsen vive GA, 0 6, 9228 om 228 em 8 GHz with 40 MHz frequency step. 
Nodes e^ enterfice hirve a ipacing of A, /20 3 M EE 







Wy end) Grection. 
121 « 122 node: on the interface. 
2 






Age @8 Ge (highart frequency of interent) 
N F 
Loop is centered N | Z 
at origin (0,0,0) 
and tilted by 45° 
about x axis 
Loop diameter: 
3 feitd cellys.625 mm) 
interface scattered fields obtained using 
modified OM approach and FEKO. ( Obs. Une 1 
Interface is parallel to yz plane at a distance of 4- corresponds to 220 
0.75 cm fromthe center of loop along x direction. Incident modulated Loop diameter: 5.625 mm. 
. f . . gaussian pulse Loop cross section radius = x of the loop circumference. 
Results along the observation line are compared at 4 GHz which is the center 
MoM results agree well with FEKO results. 


frequency of the modulated Gaussian pulse. 


MoM code vs FEKO results comparison E 
Observation line 3 


Obs. Line 3: Interface E, field comparison Obs. Line 3: Interface E, field comparison 





Magnitude Phase 











Line 3: x = 0.75 cm, Z = -7.5 cm 





Tilted Circular Loop 


In the frequency domain, the modulated 
Gaussian pulse has a range of 40 MHz to 
8 GHz with 40 MHz frequency step. 


Magnitude Phase 



































trnertece sine Ch 0 6, 2228 om 0 228 em 
Nodes e^ enterfice hirve a ipacing of A, /20 3 Plone weve . 
inymi: drenien. Theta = 45, phi=45 deg. 
121 « 122 node: on the interface. 
———————Ó 2 
” Fa 
P 
Loop is centered p 
at origin (0,0,0) 
and tilted by 45° 
about x axis 
Loop diameter: 
3 fdtd celis5.625 mm) 2 
interface scattered fields obtained using 
modified DM approach and FEKO. ( Obs. Une 1 
we wo interface is parallel to yz plane at a distance of 4- corresponds to z=0 
0.75 cm from the center of loop along x direction. Incident modulated Loop diameter: 5.625 mm. 
gaussian pulse Loop cross section radius © x of the loop circumference. 


Results along the observation line are compared at 4 GHz which is the center 


frequency of the modulated Gaussian pulse. MoM results agree well with FEKO results. 
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Multi-scale problem: ‘Transmit Case 


The modulated Gaussian pulse spectrum has a 
range of 40 MHz to 8 GHz with 40 MHz 
frequency step 


Finite pec plate 
(21.75 cm x 
21.75 cm) 


Tilted loop antenna 
7 Obs. Line 


y 


X 


(-1.5 cm, O0|cm, 0 cm) 





Modulated Gaussian Y = at 
pulse spectrum Delta gap 


FDTD 
region 


Loop is centered at origin, and MoM 
has a diameter of 5.625 mm region 
Tilt angle of loop is 45° Interface 12cm 





Loop cross section radius = — of the loop circumference (22 .9 cm x 


Interface has 121 x121 nodes 22.5 cm) 
Nodes on interface have a spacing of 

Anig/ 20 — 0.1875 cm 

Angy € 8 GHz (highest frequency of interest) 

Ànign 7 3-75 cm 

scattered fields are computed in the time-domain along the 


observation line, and transformed to the frequency domain. 





Hybrid FDTD vs FEKO results 
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The modulated Gaussian pulse spectrum has a 
range of 40 MHz to 8 GHz with 40 MHz 
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Ez phase (degree) 
8 8 


“200 





“1 B 0.6 DA 
Distance along obs. line (X) 


Finite pec plate 


frequency step (21.75 cm x 
Tilted loop antenna 21.75 cm) 
Z Obs. Line 
a X faa) cm, 0 cm) 
i modulated y we 
E Gaussian pulse spectrum Delta gap i- FDTD 
x | region 
u Loop is centered at origin, and MoM 
= has a diameter of 5.625 mm j 
region 
Tilt angle of loop is 45° 8 Interface 12 cm 
zu 14 = s Da D4 MO ET PT m n" 24 = a Loop cross section radius = = of the loop circumference (22.5 cm X 
Distance along obs. hne (X) Distance along obs. line (1) Interface has 121 x121 nodes 22.5 cm) 


Nodes on interface have a spacing of 
Ahign’20 = 0.1875 cm 

Anion @ 8 GHz (highest frequency of interest) 
Anish = 3.75 cm 


0. Number of operations using LET. method at interface Š 
Efficiency factor = = ———_1__—W___1______ > 10 


Number of operations using hybrid method 


Hybrid method and FEKO results have good comparison. 
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Multi-scale problem: Receive Case 


In the frequency domain, the modulated 
Gaussian pulse has a range of 40 MHz to 
8 GHz with 40 MHz frequency step 


Tilted loop antenna 


? al 
Incident field: X ; > en 
0 2 90° 

p = 0° i 
Modulated Gaussian pulse y 


Loop is centered at origin, and MoM 
has a diameter of 5.625 mm : 
Tilt angle of loop is 45° region 


1 
Loop cross section radius = Un of the loop circumference 


Interface has 121 x121 nodes 

Nodes on interface have a spacing of 
Anig/ 20 — 0.1875 cm 

Angn € 8 GHz (highest frequency of interest) 
Arion = 3.75 cm 


hig 





Finite pec plate 
(21.75 cm x 
21.75 cm) 


Obs. Line 


: Dem) 


FDTD 
region 





Interface 
(22.5 cm x 
22.5 cm) 


Scattered fields are computed in the time-domain along the 
observation line, and transformed to the frequency domain. 


Ez magnitude (V/m) 
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Hybrid FDTD vs FEKO results 
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0. Number of operations using LET. method at interface Š 
Efficiency factor = ————————————— > 10 


Number of operations using hybrid method 





4 GHz 


900 


a ^ 400 
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-1 -0.75 -0.5 
Distance along obs. line (A) 


-0.25 -0.8 -0.6 -0.4 -0.2 
Distance along obs. line (A) 


In the frequency domain, the modulated 
Gaussian pulse has a range of 40 MHz to 
8 GHz with 40 MHz frequency step 


Finite pec plate 
(21.75 cm x 


Tilted loop antenna 21.75 cm) 


! 
Z -Al 
i | n 
y e 


MoM 
n Interface 
Loop cross section radius = —- of the loop circumference 


Interface has 121 x121 s (22.5 cm x 
ace nas X nodes 
22.5 cm) 


Nodes on interface have a spacing of 
Agigy/20 — 0.1875 cm 

Ag, (0 8 GHz (highest frequency of interest) 
Anish = 3.75 cm 


Incident field: 

9 — 90° 

®=0° 

Modulated Gaussian pulse 


om, 0 cm) 


FDTD 
region 


Loop is centered at origin, and 
has a diameter of 5.625 mm 
Tilt angle of loop is 45° 





12cm 


Hybrid method and FEKO results have good comparison. 





PENNSTATE 





Handling arbitrary shaped surfaces 























Up 
X3,Y3,23 " Observation point. ^ — Observation point. 
Apex 2 Consider a k p j P 
- currentpath  / K>,Y>,Z5 X A j 
= / Common X3,Y3,23 , 4 R pam » N9 Apex 2 / f 
Common edge 2 | edge point Apex 2 Na. n D A Y, | ff 
aa XUYi4 ^, / / 
* / t/ J 
Current path Z Apex 1 ye a / i lw 
O F XYL P i d f 
(origin) - Apex 1 ; ~Moditied RWG A 
Modified RWG basis function 8 / 4 CB wv 
"EP 
Example: Sphere discretized basis function lacur pan P 
: Modified RWG 
into triangular patches Apex i sci cedi ier 1 
current path. 
(a) Current distribution (b) Current distribution 
y s : à in top part of current in bottom part of current 
Scattered field expressions from basis function current  pathon modified ewe path on modified RWG 
basis and distance basis and distance 
Ona current path: parameters for an parameters for an 
observation point. observation point. 
e jBRi -jfr z 1 1 
Eam = Uni n= cos(BH ) — ju,sin(BH )e /"" (| +- . 
ut = —J30Im | p OPE Geena ee as I= Imsinfp* 




















Ir. "uet e Jr jsin(BH ) 
n = (u A wycos(pn-)— PT ripe pul +o?) | | 
d pr Bent current path is handled using a procedure 
-jBR> „-—jßr similar to the one used for bent wires to obtain 
e 2 e — 1 
Ey. = -j PR m cos(BH*) ju; sin(BH*)e JP (7 j a3) time compatible scattered fields. 
M Brant NEUEM ) S 
Ey2= (uz + aa s e JP" (rBu$ 4 jvj 
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NUMERICAL RESULTS 


Hybrid FDTD method for problems involving 
wire and surface scatterers 
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Problem definition 





Validation of the MoM code 


Geometry is centered at origin (0,0,0) 
Plate side length = A/2 
Wire length = 3.75 mm (2 FDTD cells) 


Interface size: GA x Aia 7 22.5 cm x 22,5 em 


Nodes on interfacehave a spacing of Au /20 
in x and y directions. 
121x121nodesontheinterface. 


Mim © 8 GHz [highest frequency of interest) 
Voltage gap is located at 


Modulated Gaussian wire center 


pulse spectrum a- P Wire 


The modulated Gaussian pulse spectrum has a | Plate 
range of 40 MHz to 8 GHz with 40 MHz 
frequency step 
Interface scattered fields obtained using 
modified RWG approach and FEKO. 





Interface is parallel to xy plane at a distance of A/10 
from the scatterer top point along *z direction. 


Voltage gap used for excitation 


1, 
Wire radius = <p Wire length 





,* A10 


Obs. Line 1 corresponds to x=0 


Double precision accuracy 
Plate discretized into 262 triangles, option chosen in FEKO 
triangle edge length =A/20 


Wire discretized in 2 segments 


Scattered fields are computed in the time-domain along the 
observation lines, and transformed to the frequency domain. 
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MoM code vs FEKO results comparison 
Observation line | 








E, field comparison along observation line 1(x=0, z=0.75 cm). E, field comparison along observation line 1(x=0, z=0.75 cm). 
Magnitude Phase Magnitude Pise 
-4 = 
—— modified RWG —— modified RWG = = =FEKO 
Nu 
E @ 100 E F iga 
> 4 o 2 p 
o o 9 0.06 S 
"UO c = 
23 g 0 t g 9 
& S 2 0.04 £ 
E 2 = = N 
x 
ui 4! -100 H 0.02 zum 
1 
- 0 -200 
7o 50 100 150 ai 50 100 150 0 50 100 150 
Interface Node Interface Node Interface Node Interface Node 
Geometry is centered at origin (0,0,0) as adit 
Interface size: 6 i X GAgig, = 22.5 cm x 22.5 om Plate side length - M2 — A em 
Wire length = 3.75 mm (2 FDTD cells) 
Nodes on interfacehavea spacing of Ain, /20 i 
in xandy directions. Wire radius = — wire length 
121 x 121 nodeson theinterface. 50 
A eer OMIT Voltage gap is located at 
Modulated Gaussian wire center 
pulse spectrum “A Wire . $ M10 
The modulated Gaussian pulse spectrum has a Plate 
range of 40 MHz to 8 GHz with 40 MHz 
MS step Obs. Line 1 corresponds to x-0 
Interface scattered fields obtained using P " 
modified RWG approach and FEKO. T—— Double precision accuracy 
PP Plate discretized into 262 triangles, option chosen in FEKO 
Interface is parallel to xy plane at a distance of A/10 triangle edge length - A/20 
from the scatterer top point along *z direction. Wire discretized in 2 segments 
Results along the observation line are compared at 4 GHz which is the center Modified RWG results agree well with 


frequency of the modulated Gaussian pulse. FEKO results. 





MoM code vs FEKO results comparison 


Observation line 2 


E, field comparison along observation line 2 (x=3.75cm, z=0.75cm). 





Magnitude Phase 
3 
jM 200 ix x10? 
- - -FEKO 
— 0.8 4 
E = 100 = 
= E E35 
Š 0.6 © 2 r 
= 9 0 2 
$04 B e 2.5 
c a © 
x M E 2 
Ww B 
02 HN Ù 
1.5 
0 -200 1 
0 50 100 150 0 50 100 150 0 50 


Interface Node Interface Node 


Geometry is centered at origin (0,0,0) 
Plate side length = A/2 
Wire length = 3.75 mm (2 FDTD cells) 


Interface size: 6 i X GAgig, = 22.5 cm x 22.5 om 


Nodes on interfacehavea spacing of Ain, /20 
in x and y directions. 


: . i 
Wire radius — — wire length 
121 x 121 nodes on theinterface. 50 


Nis € 8 GHz [highest frequency of interest) 


Modulated Gaussian 
pulse spectrum A Wire 


The modulated Gaussian pulse spectrum has a 
range of 40 MHz to 8 GHz with 40 MHz 


frequency step 
Interface scattered fields obtained using 
modified RWG approach and FEKO. 





Voltage gap is located at 
wire center 


Plate 


Obs. Line 1 corresponds to x=0 


Plate discretized into 262 triangles, 
triangle edge length = A/20 
Wire discretized in 2 segments 


Interface is parallel to xy plane at a distance of A/10 
from the scatterer top point along +z direction. 


Results along the observation line are compared at 4 GHz which is the center 
frequency of the modulated Gaussian pulse. 





Interface Node 
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E, field comparison along observation line 2 (x=3.75cm, z=0.75cm). 


Magnitude Phase 


200 


—— modified RWG 
=- =- -FEKO 


Ez phase (degree) 
o S 


8 





0 50 100 150 
Interface Node 


Voltage gap used for excitation 


$ X10 


Double precision accuracy 
option chosen in FEKO 


Modified RWG results agree well with 
FEKO results. 








MoM code vs FEKO results comparison 


Observation line 3 
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E, field comparison along observation line 3 (x=7.5cm, z=0.75cm). — E; field comparison along observation line 3 (x-7.5cm, z-0.75cm). 


Magnitude Phase 


Magnitude Phase 


200 x10? 


N 
A 


—— modified RWG 
- - -FEKO 


= 

S 
N 
N 








200 


—— modified RWG 
= = = FEKO 


—— modified RWG 
- - -FEKO 





Š 8 € g 1 
— - — o 
o P 2 2 > 
3 = O o 
Uo — 
5 g ° 218 g 0 
© g = © 
E a m = 
x x E 1.6 
ui w -100 N Ù 400 
1.4 
-200 
i aee PE m 2 = lim 1 2 50 100 1900 50 100 150 
. Interface Node Interface Node Interface Node 
Geometry is centered at origin (0,0,0) "T - 
i - oltage gap used for excitation 
Interface size: 6. x 6i, 7 22.5 cm x 22.5 cm Plate side length = M2 ge gap 


Wire length = 3.75 mm (2 FDTD cells) 
Nodes on interfacehave a spacing of i /20 


in x and y directions. 


i : i 
Wire radius — — wire length 
121 x 121 nodes on theinterface. 50 





Nan @ 8 GHz (highest frequency of interes) . 
Voltage gap is located at 
Modulated Gaussian wire center 


pulse spectrum d Wire 





The modulated Gaussian pulse spectrum has a 
range of 40 MHz to 8 GHz with 40 MHz 
frequency step 
Interface scattered fields obtained using 

modified RWG approach and FEKO. 


Plate 


Obs. Line 1 corresponds to x=0 


Plate discretized into 262 triangles, 


Interface is parallel to xy plane at a distance of A/10 triangle edge length =/20 
from the scatterer top point along +z direction. Wire discretized in 2 segments 


Results along the observation line are compared at 4 GHz which is the center 
frequency of the modulated Gaussian pulse. 


,* M10 


Double precision accuracy 
option chosen in FEKO 


Modified RWG results agree well with 
FEKO results. 
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Multi-scale problem 
Helical antenna with ground plate 


Examples with human model 









External 
Internal antena 
antenna 
PEC plate 
21.75cmx 21.75cm FDTD region MoM region 
























Internal 
antenna 
. x: (Cross-dipole) 
region 12 cm j | 








o P f 
: i . 
o © M 
Y 
" X \ 
Ss * 
IN x 
k 
A \ y% 
Ra ^ X 
à E Ñ 
L-—— X 
> i " 
x 
` 


E External b> 
a nte nn a | a á 


di "Truncated 
| torso 


- Interface 
22.5cm 


'""^—————  . The helix has 4 turns, a base diameter of 
3.62 cm, and a height of 1.2 cm. 


1.2 cm 


MoM region 

+ The cross section radius of helical wire 1s 
1 . . 

y Helix with ground plate Tog Of the helix height. 


Excitation: delta gap source with modulated 
Gaussian pulse spectrum (40 MHz —8 GHz). 7 


X Ground plate diameter = 10.2 cm. 


Scattered fields are computed in the time-domain along the 
observation line, and transformed to the frequency domain. 
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Hybrid FDTD vs FEKO results at 1 GHz 


2,5 
= 2 
E i 
= k 
o 1.5] D 
s a 
= g 
m4 = 
E = -150 
A Ww 
ui 


e 
in 






ü i = 
0.05 O06 OOF 0.08 O09 01 011 
Position (z) along obs. line 


-200 - 
0.05 0.06 2007 008 009 01 011 
Position (z) along obs. line 


E, comparison 


0.04 r- —€——— - 250 


— FEKO 
- - - Hybrid FOTO 


2 
2 


E g 200 
> 2 
v T 
3 z 
= 0.02 g 150 
S s 
s E 
i 0.01 iJ 100 


50 — 
0.05 006 007 008 009 01 011 
Position (z) along obs. line 


D = 
0.05 0.06 O07 0.08 O09 01 011 
Position (z) along obs. line 


E. comparison 


E Number of operations using I.F. T. method at interface 5 
Efficiency factor = ————— —— ——————— —————2— ——7—— 10 


Number of operations using hybrid method 






—FEKO 
- - - Hybrid FDTD 





Ey magnitude (V/m) 
S 


Ey phase (degree) 





-200 7 
0.05 006 O07 0.08 009 OFF 02011 
Position (z) along obs. line 


Ü = d i = 
0.05 006 OOF O08 O09 01 011 
Position (z) along obs. line 


E, comparison 


PEC plate 
21.75cmx21.75cm 





12 cm 


—— fobs 


"m Interface 


22.5cm 
1.2cm «———— .— "Thehelix has 4 turns, a base diameter of 
3.62 cm, and a height of 1.2 cm. 
MoM region 
Excitation: delta gap source with modulated t The cross section radius ofhelical wire is 


j z 
Gaussian pulse spectrum (40 MHz — 8 GHz). bey Helix with ground plate = of the helix height. 


X Ground plate diameter = 10.2 cm. 


Hybrid method and FEKO results have good comparison. 
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Hybrid FDTD vs FEKO results at 4 GHz 






| Eno | 
1 = = = Hybrid FOTO 
t T 
= E oj 
o F 
3 2 
5 2 -100| 
m = 
= a 
wc M 
Lu 
m -200 








0.06 007 008 009 01 011 012 
Position (z) along obs. line 


Oo 1 4 L 5 E A 4 
006 007 008 O009 01 O11 012 
Position (z) along abs. line 


E, comparison 


——FEKO 
- - - Hybrid FDTD 






0.1 


Ez magnitude (V/m] 


o 
a 





0.06 0.08 0.1 0.12 | 0.06 0.08 
Position (z) along obs. line 


0.1 0.12 
Position (z) along obs. line 


E. comparison 


Number of operations using I.F T. method at interface 


Effici factor = 
RED Number of operations using hybrid method 


10? 





0.012 100 s : ; 
—FEKO 
0.01 0 ia 
E T 
> 0.008 2 400 
i F 
& 0.006 g -200 
S 8 
E 0.004 S -300 
=h Ww 
Lu 


0.002 -400 





Ü 4 L 4 i E 4 L 
0.06 O07 0.08 009 01 011 0.12 
Position (z) along obs. line 


006 007 008 009 01 011 012 
Position (z) along obs. line 


E, comparison 





PEC plate 
21.75cmx21.75cm 
12 cm 
cm Interface 
22.5cm 
1.2cm «———— .— "Thehelix has 4 turns, a base diameter of 
3.62 cm, and a height of 1.2 cm. 
MoM region 
Excitation: delta gap source with modulated t The cross section radius ofhelical wire is 
Gaussian pulse spectrum (40 MHz — 8 GHz). E 


1 . P 
Uy Helix with ground plate Too Of the helix height. 


X Ground plate diameter = 10.2 cm. 


Hybrid method and FEKO results have good comparison. 


Ex magnitude (Vim) 
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Hybrid FDTD vs FEKO results at 8 GHz 


——FEKO 
= = = Hybrid FDTD 






Ex phase (degree) 
db Ł i 
s 8 


& 
: 





0.06 007 008 009 01 O17 O12 


Position (z) along obs. line 


006 007 0.08 009 0.1 0.11 O12 
Position (z) along abs. line 


E, comparison 









0.035 ; z r - 
= | 
0.03 > = = = Hybrid FDTD = = = Hybrid FDTD 
= SUf 
E 0.025 T 
> 2 
3 0.02 - = -100 | 
= g 
2 E 50 
N 0.01 " 
-200 f 
0.005 
0 — -250 —— — - 
0.06 0.08 0.1 0.12 0.06 0.08 0.1 0.12 
Position (z) along obs. line Position (z) along obs. line 
E, comparison 
E Number of operations using I.F. T. method at interface 5 
Efficiency factor = ————— —— ——————— —————2— ——7—— 10 


Number of operations using hybrid method 


Ey magnitude (V/m) 
= 2 = 
E & S 


2 
S 





PEC plate 
21.75cmx21.75cm 


——FEKO 
= = = Hybrid FDTD 








gy 200 
o 
a 
= 400 
d 
= 
cL 
= 
W -600 
0.08 007 008 009 01 011 012 006 007 008 009 01 01! 012 
Position (z) along obs. line Position (z) along obs. line 
E, comparison 


12 cm 





"m Interface 


22.5cm 
1.2cm «———— .— "Thehelix has 4 turns, a base diameter of 
3.62 cm, and a height of 1.2 cm. 
MoM region 
Excitation: delta gap source with modulated t The cross section radius ofhelical wire is 


Gaussian pulse spectrum (40 MHz — 8 GHz). E 


1 . P 
Uy Helix with ground plate Too Of the helix height. 


X Ground plate diameter = 10.2 cm. 


Hybrid method and FEKO results have good comparison. 





NOVEL FINITE DIFFERENCE METHOD 


Capacitance matrix computation for coated 
multi-conductor problems 
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Proposed Finite Difference method 


2D stencil for finite difference method 


: Laplace’s Equation: V*@ = 0 
e 
à ly 2) 4- 0(3) 4 9(4) 4 9(5 
Pla a t sy -22 OK Ø(5) 
4 ÂX 
9 


Handling non — cartesian geometries: 


* The potential at node 2 is expressed in terms of potential at nodes 1, 3, and 4 (refer figure) 


Proposed boundary condition 
Constant K1 Constant K2 Constant K3 


The computational 
domain is truncated as 


| 7 
TOE (a) (1) + (| (3) - (Goce (4) to fo shown, let us consider 











(Ax Ap) (2Ax Ap) (AxAp) (Ax) (24x (Ax) ie er n T is 
fictitious. 
* An = 0,1 here when considering node 2. At the boundary node 1: 


* pija PE) OE) GG) OG) 


We assume a fictitious node 5 
and express the potential at 
fictitious node 5 in terms of 
potentials at nodes 1, 4, and 6 
using quadratic extrapolation. 
Therefore, we have at node 1: 


* Therefore, (2) can be expressed as: 


* $(2) 2 k19(1) - 293) 4 k34(4) 


All the nodes in this equation are é(1)- (2) +(3)+(4) +19 (1) +e2(4)+e3p(6) 
inside the computational domain 4 





* We express the potential at all such special nodes in the 2D FD grid in a similar manner. 


PENNSTATE 





Handling dielectric coated structures 


Step 1: We express the potential at node 1 in terms of the potential at nodes 1*, and the 
known potential ‘V’. Boundary condition at node 1* is used. 


Dielectric coating 
with thickness t 
Conductor 


8 cells Ap: distance of node 1 from 1* 


Ax: FD grid spacing between 
nodes 





Ground: OV 
V: potential of conductor 


g Relative dielectric constant 
of the dielectric coating 


An illustration of dielectric coated structure 


A staircase approximation of the 
dielectric coated metal structure is 
assumed as shown in the picture 


*. D; : displacement current vector in free space at node 1" 

* D,:displacement current vector in dielectric at node 1* 

*. d(1"): potential at node 1*: (1): potential at node 1; V: potential of the metal/conductor 
* Boundary condition at node 1" : 





— E E E field unit vectors at 
— | a | ary node 1* ! an oat 
* (D; -D, ). fi 20 PORE , boundary node 1* in free 
| J space and dielectric resp. 
. COUR LU). 
(Ap) 
* Rearranging equation (2), we get: Constant C. 
Constant C, i 
. (1° = [c y a. ( a 
p(1") = (<p) OD t (zov) 


* So, we have: b(1°) = Gol) + G 
Step 2: We express the potential at node 1* in terms of potential at nodes 1, 2, 
and 3 using quadratic interpolation. 
* The potential at node 1 is expressed in terms of potential at nodes 1 a2 2, and 3. 
Constant il Constant K2 


/ 
(Axy(24x) 


/ 
| . (dp)(24x) (Ap)(Ax) 
' p(t) epl )* aras) 9 )- (set (3) 
* Therefore, d( 1) can be expressed as: 
* $(1) 2 kip(1*) * k2(2) *- k39(3) 


* C1) = RCC, PCL) + Ca) + K2p(2) + k3p(3) 


* We express the potential at all such ‘special nodes in the 2D FD grid (nodes 
adjacent to dielectric interface)’ in a similar manner. 
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Handling two layer dielectric coated structures 


(0 c-2 a M *—3 FD dormain grid lines 


Dielectric 2 (€,,) 





Dielectric 1 (€,,) 
Metal at known potential 


O(P)-9(Q) _ : 0(Q)—O(R) 
DU =“ 


rl t1 r t2 
ØCQ)-ØCR) | 0(R)-9(1) 
r2 t2 Ap 


Ø(R)=aØ(1) +b 


E (4x) (24x) (4p) (24x) (Ap)(Ax) 
HO = (arr apea ca) 9 * (ce amaro) 9 - [a5 9 





(Ap)(Ax) 


E (4x) (24x) (4p)(24x) 
90 - (naa c2) OO +D + (aer aaa) O (C0) tA 
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Three trapezoid shaped conductors with two layer 
dielectric coating: Problem definition 


60 cells 





An example of a practical problem 
in packaging applications. 





tic coating 1 (er = 5) 





42 cells 





Narrow channel 


> 


12 cells 





Dielectric coating thickness = 0.1 units 








OV 


Handling narrow channel region 


V¥=0 boundary condition 1s used | 
Dielectric coating 2 
Dielectric coating 1 : i 


on all four sides 


Coated parallel plate capacitor 





Metal 





Narrow channel region represented as 
staircased coated parallel plate system. 


Potential and Electric Field distribution 
LV, OV, OV excitation of conductors 


Potential distribution 


ages Potential distribution, V(x,y) in volts 







y-axis in meters 
o 





-20 -10 0 10 20 30 
x-axis in meters 


Contour plot 





Dielectfic coating 1 (er - 5) 
42 cells 
Dielectfic coating 2 (er = 9) 


Dielectric coating thickness = 0.1 units 


V—0 boundary conditionis used 
on all four sides 





1: 
3r 
E Contours 


al | | | | | L1 | | i} Í | 
:292821-26252423222*2019181716151413121110.9 .8 .7 48:5 4.3.2-10 1224 5 8 7 8 910111213141516 1718192021222324252827 282930 
celis along x 
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Electric field distribution 





0.18 


0.14 


0.12 











Cells along x —> 


Electric field at nodes on contour enclosing the 
conductor is used to find enclosed charge. 


2 D. fi = Gene. 


C1101 * C150; t C450, 
C5401 + C2202 + C2303 
C3101 + C3202 + C3303 


Qı = 
Q2 = 
Q3 = 


Potential and Electric Field distribution 
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OV, IV, OV excitation of conductors 


Potential distribution 


arene Potential distribution, V(x,y) in volts 


0.8 

uv 

5 

3 0.6 

E 

£ 

o 

Š 0.4 
0.2 





-20 -10 0 
x-axis in meters 


10 20 30 


Contour plot 





60 cells 






8 cells 
«—— ———- 


Dielectfic coating 1 (er - 5) 
42 cells 


Dielectfic coating 2 (er = 9) 


Dielectric coating thickness = 0.1 units 


V—0 boundary conditionis used 
on all four sides 


Electric field distribution 


iem h m Hour m orm 







LLE Ed 


aulcm mau ELA Rim michnhe m 


HTTP RO dH 


p i 


Cells 
alon 


j =i 


Y 





b 5i Contou rs a 


11.1 
nn eee bee dd Kha tad ee Po 


ls ag 


Cells along x | —« 


Electric field at nodes on contour enclosing the 
conductor is used to find enclosed charge. 


» D. fi = enr, 


Q1 = C1101 + C1202 + C1303 
Q2 = C2101 + C2202 + C2303 
Q3 = C310, + C3202 + C3303 
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Capacitance matrix results 


Resu Its Com Da rison All C values are in pF 


2 CNN RN 





FD (coarse mesh + interpolation) 48.0 59.35 48.0 -22.97 -22.97 -0.539 

FD {2x fine mesh + interpolation) 48.5 59.95 48.5 -23.2 -23.2 -0.55 
Ansys 2D 49.06 61 49.06 -23.8 -23.8 -0.572 
Ansys 3D 50.37 63.07 50.77 -24.87 -24.83 -0.2084 


C33 = C, C = C35, C5 = C4 using symmetry 
FD results obtained using coarse mesh 


60 cells . i 
ERN and interpolation have good accuracy. 


OV 


Maximum percentage difference 
ic coating 1 (er = 5) between FD results and Ansys 2D 
ic coating 2 {er = 9) results = 2.7%. 






42 cells 


Dielectric coating thickness =. 1 units 





Ansys 3D results are converged 
results, for a length of 3m. 


V) boundary condition is used 


on all four sides > . 
Time Resource comparison 
E a E A T pe oarsemesh | axfine mesh 
Efficiency factor = Time 10810 sec. 44280 sec. 


ee 
Time resource utilized using coarse mesh and interpolation 
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Observations and Conclusions 


Magnitude and phase of the scattered field obtained using the proposed hybrid method for multi- 
scale problems demonstrated are in close agreement to that obtained using commercial method of 
moments code. 


The hybrid method provides accurate results for the transmit as well as receive mode of operation. 


Hybrid results have best accuracy at the center frequency of the band (B = B,) , and a range for 
which the accuracy is found to deviate by about 10% from the accuracy at center frequency is 
found. 


Results have good accuracy when scatterer (arbitrary wire or wire and surface combination 
problem) spans multiple FDTD cells. 


The proposed hybrid method is much more efficient (factor>10*) when compared to the past 
approach where Inverse Fourier Transform (I.F.T.) 1s used at the interface to obtain time domain 
fields used as source for the FDTD region. 


The proposed finite difference algorithm provides good accuracy and efficiency (factor >4 when 
compared to 2x fine mesh) for capacitance matrix calculation for dielectric coated, multi-conductor 
problems where conductor boundaries do not align with grid lines. 
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Future Research 


* Extend the hybrid method to utilize the near field region for arbitrary shaped wires 
and surfaces. 


* Extend the finite difference code to 3 dimensional problems when longitudinal 
dimensions are finite. 
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